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Stellar structure in braneworlds is marlcedly different from that in ordinary general relativity. As 
an indispensable first step towards a more general analysis, we completely solve the "on brane" 
4-dimensional Gauss and Codazzi equations for an arbitrary static spherically symmetric star in a 
Randall-Sundrum type II braneworld. We then indicate how this on-brane boundary data should be 
propagated into the bulk in order to determine the full 5-dimensional spacetime geometry. Finally, 
we demonstrate how this procedure can be generalized to solid objects such as planets. 
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I. INTRODUCTION 

One physically important problem in the braneworld 
scenario is the development of a full understanding of 
stellar structure and black holes |1„ 2, 3, 4, 5j. What 
is already known is that stellar structure in braneworlds 
is rather different from that in ordinary general relativ- 
ity. The key point is that if one is confined to making 
physical measurements on the brane, then the restricted 
"on-brane" version of the Einstein equations does not 
form a complete system for specifying the brane geome- 
try V, "s] . The problem lies in the fact that the bulk 5- 
dimensional Wcyl tensor feeds into the "on-brane" equa- 
tions and so connects the brane to the bulk. 

On the other hand, one may solve the "on-brane" 
version of the Einstein equations to obtain a consistent 
class of boundary data that satisfies the 4-dimensional 
Gauss and Codazzi equations ^ , analogously to the ini- 
tial data in the standard (3 + 1) decomposition of glob- 
ally hyperbolic manifolds in general relativity. Although 
the relevant partial differential equations are now elliptic 
rather than hyperbolic, one can use the boundary data 
in a completely analogous fashion — as input into the 
5-dimensional bulk Einstein equations in order to "prop- 
agate" the 4-dimensional geometry off the brane and into 
the bulk [113. 

In this Letter we present an algorithm for completely 
solving the 4-dimensional Gauss and Codazzi equations 
for a static spherically symmetric star on the brane. This 
provides the most general boundary data, suitable for 
then determining the bulk geometry by an appropriate 
relaxation method 0. 

Specifically, we will use Gaussian normal coordinates 
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adapted to a timelike 4-dimensional hypersurface (space- 
like normal) in a 5-dimensional geometry. 



ds5 = drf + gab da;" da;^ 



(1) 



and assume a type II Randall-Sundrum braneworld 0, 
in which the bulk metric is a 5-dimensional Einstein 
space, (i.e., with Ricci tensor proportional to the met- 
ric), but the 4-dimensional Lorentzian signature metric 
gab is as yet undetermined. 

If we impose ^2 symmetry on this spacetime, and tune 
the 5-dimensional bulk cosmological constant and brane 
tension appropriately, then the junction conditions to- 
gether with the projected 5-dimensional Einstein equa- 
tions (the Gauss equations) reduce to the induced but in- 
complete on-brane "Einstein equations" in the form 0,0 
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where is a constant inversely proportional to the brane 
tension, A4 is an on-brane cosmological term, the nonlo- 
cal term Sab is simply the projection of the 5-dimensional 
Weyl tensor onto the brane, and in addition to the 
usual stress-energy term there is also a nonlinear term 
quadratic in stress-energy. In addition to ^ one finds 
that the extrinsic curvature 



1 dgab 

2 dr] ' 



is related to the on-brane fields via 
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The Codazzi equation 

WbiK'''' ^ Kg"''] =0, 
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is then equivalent to 4-dimensional stress-energy conser- 
vation for the on-brane matter: VTab = 0. Eqs. Q 
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and Q complete the specification of the boundary data 
by effectively supplying the "on brane" metric and its 
normal derivative. 

The only truly general thing we know about the nonlo- 
cal Weyl tensor projection term £ab is that it is traceless, 
f^a = 0. Together with |5J) this implies that 

R = -Svr T + 4A4 - [T -T- \T^] . (6) 

In the vacuum region outside a star this reduces to 

i? = 4A4, (7) 

whereas inside the star it takes the form 



R = [nonlinear-source] . 



(8) 



If for instance we are dealing with a perfect fluid then 

R = 87:{p - 3p) + 4A4 + - 3p2) + i(p - 3p)^] . 

(9) 

If we are restricted to making physical measurements 
on the brane, then ^ is the only really general thing 
we can say. Eqs. (0), © are much weaker than the 4- 
dimensional Einstein equations and so the solution space 
will be much more general. Of course, we must also 
use the extrinsic curvature constraint Q to complete 
the boundary data, and then ultimately probe the bulk 
geometry off the brane, and this will indirectly provide 
further restrictions. 

Our aim here is to solve Eqs. jSJ and in full gener- 
ality for arbitrary static spherically symmetric solutions 
on the brane, to provide appropriate data to "propagate" 
into the bulk (via Eqs. (|32|l below), to thereby test the 
consistency of candidates for realistic stellar models. 



II. VACUUM [A4 = 0] 

In this section we will consider the case when A4 = 0. 
For a braneworld star, in the vacuum region outside the 
surface we are interested in solving R = 0. In the usual 
way we can locally adopt on-brane coordinates such that 



ds^ = - exp[-20(r)] dt^ 



W) 



(10) 



If we now impose the condition i? = we have one differ- 
ential constraint connecting two unknowns — therefore 
there will be a nondenumerable infinity of solutions pa- 
rameterized by some arbitrary function of r. Various 
specific cases have already been discussed in the litera- 
ture, but no attempt has been made at extracting the 
general solution. 



A. Special cases 

Known special case solutions include: 
— "Reissner-Nordstrom-like" 



B{r) = exp[-2(/)(r)] = 1 



2M 



Z 



(11) 



Note that the parameter Z is not an electric charge, but 
should be thought of as a tidal distortion parameter. 
— "Spatial Schwarzschild" 



B{r) = 1 



2Af(l -f- e) 



exp[-2(/.(r)] = 



1 + e 



v 



(12) 



(13) 



(This geometry is also discussed in a rather different con- 
text in ref. (Iflj.) 
— "Temporal Schwarzschild" 



^ri-^^ (i 



exp[-20(r)] = 1 . 

r 



(14) 



(15) 



However these are all very specific special cases; and are 
in no way general. 



B. General vacuum solution 

Let us first write the metric in the form 

dr^ 



ds = — exp 



-2 



g{r) df 



dt' 



B{r) 



dn^ 



(16) 

The function g{r) is interpreted as the locally- measured 
acceleration due to gravity; it is positive for a inward 
acceleration. 

Now calculate the Ricci scalar R 



R = 



(2r + r'^g)B' + {2r'^g^ + 2r'^g' + Arg + 2)5-2 



(17) 

For the vacuum case we set this equal to zero, which 
yields 

(2r + r'^g)B' + {2r^g^ + 2r'^ g' + Arg + 2)B-2 = Q. (18) 

If we treat g{r) as input and view this as a differential 
equation for B{r), it is a first-order linear ODE, and 
hence explicitly solvable. The integrating factor is 

^ „^ f f ' + ^^W^y);;^V(q „,| _ 

Uro r{l+rg{r)/2) J 

(19) 

where tq is any convenient reference point, and the gen- 
eral solution is 

(20) 
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Whereas rg is an arbitrary gauge parameter, the con- 
stant B^tq) is related to physical parameters such as the 
mass and post-Newtonian corrections. In the case of the 
"temporal Schwarzschild" solution H14|l . (|15|l . for exam- 
ple, B{ro) = 2{a + ro) - 3M. 

Alternatively one may use B(r) as input to generate a 
first-order ODE quadratic in g{r) — a Riccati equation. 
However, this has no comparable general solution. The 
algorithm above is a modification of Lake's construction 
for generating spherically symmetric perfect fluid space- 
times llj. We suspect that an isotropic coordinate ver- 
sion, based on 0, may also be viable. 

This general solution can be somewhat simplified by 
using integration by parts on the integrating factor to 
remove the derivatives of g{r). We find 



Fir;ro) 



il + rg{r)/2f 
(1 + ro + giro)/2f 

,,J r ^ + rg{f)+f-g{fY 



(21) 



^ r [1 + \rg{r)\ 
ro [l + ^rogiro)]' 

where we have defined 



F2{r;ro) = exp 



Eq. H2U|) then becomes 



g{f)[l + 2fg{f)] 
2(1 -t- irg(f)) 



(22) 



(23) 



B{r) = <j / (l + f5(f)/2)F2(f;ro) df 

-ro [l + ^rogiro)f i?(ro)| 
c {r [l + irg(r)]' F2(r;ro)} 



(24) 



This is now the explicit general (static spherically sym- 
metric) solution of the equation R = using the function 
g{r) and the constant Bijo) as arbitrary input. The con- 
struction is completely algorithmic. 



J 



There is a potential subtlety if the Schwarzschild co- 
ordinate r does not increase monotonically with respect 
to the outward radial proper distance, tr- For ordinary 
stars in general relativity the monotonicity of r with re- 
spect to is guaranteed by the null energy condition. In 
braneworld stars there is no particular reason to believe 
in the monotonicity of r{£r), but our construction will 
still hold piecewise on monotonic intervals. 



III. STELLAR INTERIOR 

We now want to solve R — S{r), with S{r) a speci- 
fied source. We can now consider arbitrary values of A4 
without additional complications by including a possibly 
non-zero A4 in S{r). Proceeding exactly as above we find 



B{r) = 



l + f^S{f)/2 
f[l + fg{f)/2] 



F{f- ro) df + B{ro) 



xFir^roY 



(25) 



with the same integrating factor F(r;ro) of Eq. H21|l . 
Equivalently 



S{f)) {l + lfg{f))F2{f) df 




-ro [l + |rog(ro)]' B(ro)j 
; |r [l + ^rg{r)]^ F{r;ro) 



(26) 



Once the source S{r) is specified, this is fully general. In 
addition one must specify an arbitrary function (?(r), and 
a single arbitrary constant B{ro), and so algorithmically 
determine the metric on the brane. 

Of course the source S{r) is somewhat restricted in 
that it is an algebraic function of the on-brane stress- 
energy tensor, which is itself restricted by 4-dimensional 
energy-momentum conservation (equivalent to the Co- 
dazzi equation (jSJ as mentioned above). For a static 
perfect fluid with spherical symmetry the stress-energy 
has the form 



rpab 



p{r) exp [+2j^g{r) dr] 








p{r) B{r) 
p{r)/r^ 







p{r)/{r^sm^9) 



(27) 



The equation of energy-momentum conservation gives 

dp 
dr 



-g{r) [p + p], 



(28) 
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which, being a Hnear first-order ODE, has the exact closed-form solution 



p(r) = exp 



g{r) dr 



X i Piro) 



g{r) p{r) exp 



g{f) df 



dr 



(29) 



The lesson now is that to find all possible conserved 
tensors T"'' one is free to specify the function p{r) and 
the number p(ro) arbitrarily, and thereby calculate p(r), 
which now yields the full tensor T"^. From T°'^{r) we 
now calculate S{r), thereby fixing the intrinsic geometry 
on the brane. By eq. Q this also automatically generates 
the most general possible candidate for the extrinsic cur- 
vature K'^^ compatible with the assumed symmetries. It 
must be emphasised that only some of these braneworld 
geometries are physically meaningful, because one now 
needs to extrapolate them off the brane to see if the 
"graviton" is still bound [l3l |. 



IV. EXTRAPOLATING OFF THE BRANE 

The full algorithm is: 

• Step 1: Pick an arbitrary density distribution p{r) 
of matter inside the star; an arbitrary function g{r); 
and a single number p(ro). Calculate p{r) and so 
evaluate the stress-energy tensor T"'' and the source 
term 



S{r) = -SttT + 4A4 ~\n'^[T-T- ^T^] 



(30) 



• Step 2: Armed with this source term S{r) and the 
previously chosen g{r), pick one additional number 
B{ro) in order to calculate the function B{r) and 
thereby generate a candidate metric gab for the on- 
brane physics. 

By calculating the 4-dimensional Einstein tensor 
for this candidate metric, one can rearrange Eq. (0) 
to calculate Sab = Cna-qb and so the find the pro- 
jection of the five-dimensional Weyl tensor on the 
brane. 

• Step 3: Using this on- brane metric, and the on- 
brane stress-energy calculated in Step 1, use Eq Q 
to calculate the extrinsic curvature Kab- 

• Step 4: Evolve the metric off the brane. Sufficiently 
near the brane one can certainly use normal coor- 



dinates and so the standard result 

dKab 



rjb 



dr} 



K K'' 



applies. Rearrange this to yield: 



d'^gab 



1 „a'6' ^9aa' dgb'b 

dr] df] 
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Tjarih • 



(31) 



(32) 



Here '^^^Rnarjb is calculable in terms of the projec- 
tion of the five-dimensional Weyl tensor and the 
bulk cosmological constant. 

This now provides a well-determined set of equations for 
extending the on-brane metric into the bulk. Whether 
or not the resulting braneworld model is actually viable 
depends on whether or not the bulk geometry is "well 
behaved" — in particular, is the graviton bound or un- 
bound 10. 

We must make the caveat that the Gaussian normal 
coordinates system is likely to break down at some 
stage as one moves away from the brane. This is not 
really a fundamental objection but more of a technical 
issue. While it is easiest to set up the "on-brane" bound- 
ary conditions using Gaussian normal coordinates, as a 
practical matter when it comes to numerically solving for 
the bulk geometry one should be prepared to dynamically 
adjust the coordinate system. 



V. SOLID PLANETS 

If we are dealing with a situation of spherical symmetry 
that does not correspond to a perfect fluid, such as a 
solid planet, the radial and transverse pressures would 
be different and we would have 



R 



87r(p - pr - 2pt) + 4A4 



(33) 



1 „2 

'4 



Pr 



2pi) + \{p-pr^2ptf 



The other major change arises in the on-brane stress- 
energy tensor, which becomes 



p{r) exp [+2 g{r) dr] 



rjiab 





Pr{r) B{r) 








Pt{r)/r 



(34) 



pt(r)/(H sin 
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Then covariant conservation gives the shghtly more comphcated equation 



ar r 



(35) 



This is still a linear first-order ODE, and has the exact closed-form solution 



Pr{r) — r exp 



'l{r) dr 



X / exp 

'ro 



7(f) df 



(2 rpt{r) - g{r) p{r) r^) dr + pr{ro) > . (36) 



The lesson now is that to find all possible conserved ten- 
sors T"'' one is free to specify the functions p{r), pt{r), 
g(r), and the number Pr{ro) arbitrarily, and thereby cal- 
culate Prir) which now yields the full stress-energy ten- 
sor T"''. This is now used to calculate S{r). Once B{ro) 
is specified the on-brane intrinsic geometry (encoded in 
B{r)) is calculable. Ipso facto, this procedure also gener- 
ates the most general possible candidate for the extrinsic 
curvature K'^'' compatible with the assumed symmetries. 
Thus spherically symmetric solid objects such as planets 
are not much more difficult to deal with than are fiuid 
objects such as stars. 

VI. CONCLUSION 

What we have done in this article is to find the most 
general set of "on-brane data" suitable for characteriz- 
ing an arbitrary static spherically symmetric braneworld 
star. This on-brane data is characterized by two arbi- 
trary functions g{r), and p{r), and two arbitrary con- 
stants B(rQ), and p(ro). This is an enormous dataset, 
much more general than the various special cases dis- 
cussed in 0,0, 13 . After suitable notational modifica- 
tions, this boundary data can be used as input into Wise- 
man's relaxation algorithm for determining the bulk ge- 
ometry T'l . Although Wiseman's algorithm T] is general, 
his numerical examples were restricted to special choices 
of boundary data. Hopefully, by characterizing the most 



general set of boundary data, we are providing the ingre- 
dients for a more systematic treatment of the numerical 
problem in future work. 

Ultimately, it will be the bulk geometry that deter- 
mines whether or not a specific set of on-brane data leads 
to a physically meaningful star jl^] . We believe the tech- 
nique we have developed here could be quite powerful 
because it is completely algorithmic. Ideally, we would 
hope that broad statements about the class of possible 
stellar structure functions S{r) could made. For exam- 
ple, by a suitable characterization of possible pathological 
singularities in the bulk geometry, it might be possible 
to rule out particular classes of structure functions. 

We conclude that braneworld stars are potentially 
much more complicated than standard general relativis- 
tic stars, and emphasise that the coupling to the bulk will 
play an important role in restricting the possible func- 
tions g{r) which enter the Tolman-Oppenheimer-Volkoff- 
like equation H28() . 
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Note added: After this work was completed a paper by 
Bronnikov and Kim appeared jl4j which also derives the 
pure vacuum solution H21() ~ (|24|l . but in the context of 
braneworld wormholes. That analysis docs not consider 
the presence of matter, nor does it consider the extrinsic 
curvature. 



[1] T. Wiseman, Phys. Rev. D 65 (2002) 124007 

"arXiv:hep-th/01 11057". 
[2] N. Deruelle, "Stars on branes: The view from the brane" , 

arXiv:gr-qc/0111065 
[3] C. Germani and R. Maartens, Phys. Rev. D 64 (2001) 

124010 arXiv:hep-th/0107011 . 
[4] R. Casadio, A. Fabbri and L. Mazzacurati, Phys. Rev. D 

65 (2002) 084040 arXiv:gr-qc/0111072 . 
[5] R. Casadio and L. Mazzacurati, "Bulk shape of brane- 
world black holes", arXiv;gr-qc/0205129 
[6] T. Shiromizu, K. Maeda and M. Sasaki, Phys. Rev. D 62 

(2000) 024012 arXiv:gr-qc/9910076 . 
[7] J. Garriga and T. Tanaka, Phys. Rev. Lett. 84 (2000) 

2778 arXiv:hep-th/9911055 . 
[8] S. B. Giddings, E. Katz and L. Randall, JHEP 0003 



(2000) 023 |arXiv:hep-th/0002091| . 

[9] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 

4690 arXiv;hep-th/9906064 . 
[10] N. Dadhich, S. Kar, S. Mukherji and M. Visser, Phys. 

Rev. D 65 (2002) 064004 arXiv:gr-qc/0109069 . 
[11] K. Lake, "All static spherically symmetric perfect fluid 

solutions of Einstein's Equations", arXiv:gr-qc/0209104 
[12] S. Rahman and M. Visser, Class. Quantum Grav. 19 

(2002) 935 arXiv:gr-qc/0103065 . 
[13] A. Chamblin, H. S. Reall, H. a. Shinkai and T. Shiromizu, 

Phys. Rev. D 63 (2001) 064015 arXiv:hep-th/0008177 . 
[14] K.A. Bronnikov and S.-W. Kim, "Po ssible wormholes in 

a brane world", arXiv:gr-qc/0212112| 



